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Abstract 

We use geometric measure theory to introduce the notion of asymp¬ 
totic cones associated with a singular subspace of a Riemannian manifold. 
This extends the classical notion of asymptotic directions usually defined 
on smooth submanifolds. We get a simple expression of these cones for 
polyhedra in as well as convergence and approximation theorems. In 
particular, if a sequence of singular spaces tends to a smooth submanifold, 
the corresponding sequence of asymptotic cones tends to the asymptotic 
cone of the smooth one for a suitable distance function. Moreover, we ap¬ 
ply these results to approximate the asymptotic lines of a smooth surface 
when the surface is approximated by a triangulation. 


Introduction 

In the past decades, there has been a growing interest in defining geometric 
invariants describing singular spaces la, m, a, 0, m, 0, m, la, 

dS], dH], dO], dg. Such invariants are generally subject to two assignments: 

1. They must fit with the classical invariants when the underlying set is a 
smooth manifold or submanifold. 

2. They must satisfy some continuity conditions. For instance, if a sequence 
of singular spaces tends (for a suitable topology) to a smooth space, then 
the invariants defined on the singular spaces also tend to the smooth ones. 

After the length, the area and the volume, the most popular smooth geo¬ 
metric invariants are (sectional, Ricci, scalar, principal) curvatures, which are 
smooth functions (or tensors) defined on a (smooth) (sub)manifold and its tan¬ 
gent space. A classical approach to extending these curvatures to singular 
spaces W of a Riemannian manifold M, is to replace functions by measures 
on M. These measures are defined by integrating invariant differential forms 
over Borel subsets. As an example, Lipschitz-Killing curvature measures for 
singular spaces of a Riemannian manifold can be defined as follows: 
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• The first step is to generalize the unit normal bundle of a smooth subman¬ 
ifold to singular spaces. This has been done for convex subsets and for 
subsets of positive reach. More generally, the theory of the normal cycle 
[21] allows us to define an integral current on a large class of singular sub¬ 
sets, called “geometric subsets” of a Riemannian manifold, generalizing 
the unit normal bundle. 

• The second step consists of defining standard differential forms on the 
tangent bundle of the ambient space. 

• Finally, integrating these differential forms on the normal cycle, builds 
invariant measures, satisfying the two assignments described above. 

We remark, however, that this building does not allow us to extract charac¬ 
teristic subsets induced by a couple (W, M) (where W is a singular space of a 
Riemannian manifold M) as asymptotic directions and principal directions. 

In the framework of geometric measure theory, more specifically using the 
normal cycle theory, we use the formalism introduced in mm and [20| to pro¬ 
pose an extension of the definition of the classical asymptotic directions to a 
large class of singular spaces W. Instead of building a new “curvature mea¬ 
sure” on any couple (W, M), we associate a map that assigns a cone (called an 
asymptotic cone) of vector fields on M to each Borel subset of M. In particu¬ 
lar, if M = E^, we can reduce the target of this map to the set of quadratic 
cones of E^. Moreover, choosing a fixed scalar r > 0, we associate a fie Id of 
cones leaving in the tangent bundle TE^ with such a couple (>V,E^). The two 
assignments described above are therefore satisfied. 

• The cones associated with a Borel subset reduced to a single point of a 
smooth surface lying in E^ are nothing but the product of the classical 
asymptotic direction(s) by the normal line of the surface at this point. 

• If W/c is a sequence of singular subsets (admitting asymptotic cones), which 
tends to a smooth submanifold M, then the sequence of asymptotic cones 
associated with Wk built over a “regular” Borel subset B of M tends to 
the cone associated with M built over the same B. 

We give some applications of this building up for surfaces in E^. 

• If a plane field is defined on a singular surface of E^, we can define the 
asymptotic directions at a point m of the surface associated with a Borel 
set (for instance a ball centered at m) as the intersection of the asymptotic 
cone defined at m, with the plane defined at m, to get two directions than 
can be called asymptotic directions. By integrating these directions, we 
get asymptotic lines depending obviously on the plane field. 

• Using a triangulated approximation of a smooth surface in E^, we 
build an approximation of its asymptotic directions at some fixed points 
by considering the intersection of the asymptotic cones of the triangulation 
with the planes spanned by the faces of the triangles. We can deduce 
an approximation of the asymptotic lines of if the triangulation is a 
“good” approximation of the smooth surface. 
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Finally, we mention that this approach can be extended to other character¬ 
istic subsets or subspaces, which could be the subject of future researches. 

The article is organized as follows: 

After a summary of the main notations pertaining to Riemannian submani¬ 
folds, the first section begins with smooth surfaces in 

Since we define a new geometric invariant in quite a general context, our 
framework is an (oriented C^) Riemannian manifold M (of finite dimension). 
After a review of the geometry of the tangent bundle of M (see m for in¬ 
stance), the second section introduces the so-called asymptotic form defined in 
the tangent bundle of M, easily derived from the fundamental form [8], [20] . 

The third section uses the theory of the normal cycle [14], [21] to intro¬ 
duce the notion of asymptotic measures defined in the tangent bundle of M, 
associated with a large class of singular subsets Q (called geometric subsets in 
ED- These measures are derived from the asymptotic form defined in Section 

We describe them explicitly in classical situations (when Q is smooth or a 
polyhedron). Then, we construct the asymptotic cones. In particular, we show 
how this construction generalizes the asymptotic directions defined on a smooth 
surface in E^, fulfilling assignment above. In particular, we give an explicit 
expression of the asymptotic cones associated with a 2-polyhedron in E^. 

The fourth section deals with a general theorem of convergence of asymptotic 
cones. This is the justification of our definition, and fulfills assignment 13 In 
particular, we give explicit assumptions in terms of the fatness of the sequence 
of polyhedra, implying the convergence of the sequence of cones for a suitable 
pseudo-metric. The last section presents two applications. First, we give a 
method to build asymptotic lines on a triangulation. Then, we give a method 
to approximate the asymptotic lines of a smooth surface approximated by a 
polyhedron. We test our method on various smooth or discrete surfaces. 

Some notation - For details on the theory of smooth Riemannian submani¬ 
folds, the reader may consult for instance, [6] or m- We use the following nota¬ 
tion. Let be an n-dimensional closed (oriented) smooth submanifold embed¬ 
ded in a smooth Wdimensional (oriented) Riemannian manifold (M^, < .,. >). 
The manifold inherits a Riemannian structure by pulling back < .,. > by 
the embedding, which we still denote by < .,. >. We denote by TW^ 

{resp. TM^ M^) the tangent bundle of {resp. M^) and 
{resp. E(M^)) the space of tangent vector fields over {resp. M^). We de¬ 
note by V {resp. V) the Levi-Civita connection on (M^, < ., . >) {resp. W^). 
We denote by T-^W^ the normal bundle of the submanifold W^. The 

second fundamental form of the submanifold is the symmetric vector-valued 
(2, 0)-tensor 

h : TW^ X TW^ 

defined as follows: 

Vx G T1F^ G h{x, y) = V^y - V^y. 

Let m be a point of . The isotropic cone 

CZ' = {xeTmW^ ■.h^{x,x) = Q} (1) 

of hm is classically called the asymptotic cone of at m. For any ^ 

the eigenvalues of < hm(-, > are the principal curvatures of at m in 

the direction frn- 
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1 The case of smooth surfaces in 


We first restrict our attention to smooth closed (oriented) surfaces embedded 
in the (oriented) Euclidean space (E^, < . >) bounding a domain B. Let 

^ be the normal vector field compatible with these orientations. The second 
fundamental form of can now be identified with the tensor < h(.,.),^ > 
taking its values in C^(bE^). We denote by Ai^,A 2 ^ the principal curvatures 
of at the point m, that is, the eigenvalues of by G the Gauss curvature 
of W^, that is, the determinant of h and by H its mean curvature, that is, its 
trace (in an orthonormal frame). In a frame of principal vectors (ei^,e 2 ^) at 
m (that is, eigenvectors of hm) the matrix of hm is 

0 \ 

0 A2^; • 


At each point m G with negative Gauss curvature, the asymptotic cone 
is the union of two lines. Integrating the corresponding vector fields gives 
rise to foliations of by the so-called asymptotic curves. By definition, at 
each point, the principal normal directions of these curves (considered as curves 
in E^) are tangent to W‘^. 

Using measure theory, the goal of this paper is to define and study analogous 
cones associated with a (regular or) singular subspace of a Riemannian manifold 
M^, lying above any Borel subset of . 

Let us begin by explaining how we define such cones over any Borel subset of 
E^ in the regular case; that is, when the subset is a (compact) domain bounded 
by a smooth surface W^. Let Tv^ 2 E^ be the tangent bundle of E^ restricted 
to W^. If X is any vector field on T^sE^, we build a signed measure <I >^2 as 
follows: For any Borel subset B of E^, we write 


= / hm{prTmW^x,prT^w^x)dm, ( 2 ) 

where prT^w‘^ denotes the orthogonal projection on TmW‘^ and dm is the 
Lebesgue measure on E^. Let us now fix B and consider the map 

(3) 

where x runs over the (huge) space of vector fields E(E^)|i 4 / 2 . This map is 
quadratic in x. If we force x to be a constant vector field, then we get a 
quadratic form (that we still denote by <I>^(5)) on E^. This quadratic form has 
generically three eigenvalues, Xi{B), X 2 {B), Xs{B), that we call the principal 
curvatures of B. The corresponding eigenvectors are called the principal vectors 
of 5, and the matrix of ^^ 2 {B) in this frame is 

0 0 \ 

0 Xi{B) 0 
VO 0 Xs{B)J 

From this construction, we also introduce the isotropic cone associated with 
CP-’^ = {xGE3;$f,(B)=0}, (4) 
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(the notation coming from the fact that we restrict our isotropic cone to constant 
- that is, parallel - vector fields in E^). We call it the asymptotic cone of B (with 
respect to D). 

To clarify that this construction is linked with the classical pointwise situa¬ 
tion, suppose that B is reduced to a point {m} G W^. If ^ is a constant vector 
field such that Pm G TmW^^ then 

$^({ 771 }) = hjyi{yjjYi^yjYi^. ( 5 ) 

If z is a constant vector field such that = Cm, then 


= ( 6 ) 

This implies that in the frame (ei^, 62 ^,Cm), the matrix of $^({m}) is 

(7) 



Consequently, the asymptotic cone is nothing but the cone spanned by 

the normal Cm and . 

This construction has some advantages: One can define the asymptotic 
cones at different scales by scaling the Borel sets (for instance, by taking balls 
of radius ^ as Borel subsets). Generically, we get three geometric invariants, 
Ai(B), A 2 ( 5 ), A 3 (5), instead of two, and a two-dimensional cone instead of the 
union of two lines. Moreover, another important advantage of replacing func¬ 
tions by measures, is that this framework can be used for a large class of singular 
spaces (for instance, polyhedra, algebraic subsets, subanalytic subsets) of any 
codimension in any Riemannian manifold, as long as we can extend the notion of 
normal space. For instance, if one replaces the smooth surface by a polyhedron 
P bounding a domain P, we will get the follo wing explicit simple expression 


approximating the cone (seelSSl and also 33 for the exact formulas) 


~ ^ «{e n B)Z{e) < X, e" >2 


eGE 


and 


Qpar,!) 


{x G E^ : ^ l{e D B)Z{e) < x,e >^= 0}, 


eGE 


(see Section]^ for the notation). This is why the theory of the normal cycle, 
extensively studied over the last decades [14], [21], will be our framework. In 
the following we describe the construction of asymptotic cones in such a large 
context. 


2 Asymptotic forms 

To be self contained, we begin with a summary of the geometry of the tangent 
bundle of an (oriented) Wdimensional smooth Riemannian manifold (M^, < 
.,. >). The reader may consult [16], [17] and [ 8 ] for details. We denote by 
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TTM^ TM^ the tangent bundle of the manifold TM^. As usual, we 

consider the exact sequence of vector bundles: 


0 —> TM^ Xmn TM^ -4 TTM^ -4 TM^ Xmn TM 


N 


( 8 ) 


where i denotes the natural injection defined by 

i{ui,U2) = ^{ui ^tU2) 


t=0 


(9) 


and 


j = 

The vertical bundle of M is the subbundle V{M^) = kerj of TTM^. 
morphism i induces an isomorphism: 


( 10 ) 

The 


i : TM^ Xmn TM^ V{M^). 

If m G and x G ^ the vertical lift of ^ G at X is the vector 

= i{x^z). The morphism J = io j is an almost tangent structure on 
(J2 = 0) and V{M^) = ker J. Let 


^ : TM^ —^ TM^ Xmn TM 


N 


( 11 ) 


be the canonical vector field defined by 6{x) = (x^x) and let C : TM^ —> 
V{M^) be the vertical vector field associated with the (global) one-parameter 
group of homotheties with positive ratio, acting on the fibers of TM^. We have 
C = i o 6. We write r] = pr 2 o i , where pr 2 denotes the projection on the 
second factor of TM^ ^mn TM^ . Since is endowed with a Riemannian 
metric < •, • > and its Levi-Civita connection, we can build the corresponding 
right splitting 

7 : TM^ Xmn TM^ —> TTM^ ( 12 ) 

of the exact sequence(satisfying jf 07 = Mt^m^x^atTM^)- Let m G and 
X G . The horizontal lift of z G T^M^ at x is the vector z^ = 7 (x, z). We 

denote = Im( 7 (x, •)), from which we construct the horizontal bundle 

H{M^) such that, for all x G 

T^TM^ = ^^(M^) 0 (13) 


We denote by V(M^) {resp. 1-L{M^)) the space of vertical {resp. horizontal) 
vector fields. We denote by h : TTM^ —> H{M^) the horizontal projection, 
and by v : TTM^ —)► V(M^) the vertical projection. We remark that h = 70 J 1 . 
The morphism 

K = r]ov: TTM^ TM^ 

is the connector associated with the Levi-Civita connection. At every point 
X G TmM^ ^ the morphism 

{diTM X K), : T,TM ^ T^M^ x T^M^ 

is an isomorphism that identifies Vx{M^) with T^M^ and Hx{M^) with 
TmM^. The bundle TTM^ 1-4^ TM^ is canonically endowed with the Sasaki 
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metric < •, • > defined by the following conditions: 


( V{M^) and H{M^) are orthogonal, 

< i is an isometry, (14) 

[ 7 is an isometry. 


If 


a:M^ ^ xM^ 

is the diagonal map defined by a{m) = (m,m), then for every x G TM^^ 
{diTM X K)~^ o da{x) = 0 x^. 


Finally, the bundle TTM^ TM^ is also endowed with an almost complex 

structure F {F^ = —Id), defined by the following conditions 



II 

(15) 


1 Fh = - J. 

Therefore F\v : V{M^) - 

iJ(M^) and F\h : H{M^) - 

V{M^) are 


isometries. In this Riemannian context, we give the following definition: 

Definition 1. 1. The vector valued (N — l)-form on TM^ defined for each 

X e n{M^) by 

= [*HodgeiFCAX)]AFX (16) 

is called the asymptotic (A^ — l)-form on TM^. 

2. The vector valued {N — l)-form on defined for each x G S(M^) by 

= [*Bodge{FCAx^)]Ax'’ (17) 

is called the asymptotic {N — l)-form on . 

In this definition, *Hodge denotes the Hodge duality on for each 

X G TM^ . (The introduction of the Hodge operator in the definition of the 
generalized second fundamental form can be found in [3] when the ambient 
space is Euclidean. We adapt it here to the general Riemannian situation. It is 
equivalent to the initial definition given in [ 8 ]). Using the identification of vector 
fields and 1-forms induced by the Riemannian structure, *Hodge(^C' A X) is a 
[N — 2)-form on TL{M^). On the other hand, FX is (identified with) a 1-form, 
null on H{M^) and acting on U(M^), and *Hodge(^C' AX) is (identified with) 
an (X — 2)-form null on V{M^) and acting on 


3 Normal cycles, asymptotic measures, asymp¬ 
totic cones 

3.1 Currents and normal cycles of singular spaces 

Let Vi{TM^) be the space of l-currents of TM^ (0 < / < 2 X); that is, the 
topological dual of the space V\TM^) of I-differential forms with compact 
support on TM^ ^ endowed with the topology of uniform convergence on any 
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compact subset, of all partial derivatives of any order. The duality bracket will 
be still denoted by < •, • > if no confusion is possible. The space ViiTM^) is 
naturally endowed with the weak topology: if is a sequence of /-currents 

of TM^ and if C is a /-current of TM^, then 

lim Ck = C e V\TM^), lim < Ck^oo >=< C,uj> . (18) 

k^oo k^oo 

An /-current is rectifiable if it is associated with a rectifiable subset (see [20] for 
details). An /-current is integral if it is rectifiable and its boundary is rectifiable. 

When it exists, the normal cycle of a (compact singular) subset W of a 
Riemannian manifold is a closed integral current N(>V) G PAr-i(TM^), 
which is Legendrian for the symplectic structure on T dual to the canonical 
one on T*M^ in the duality defined by the metric. The normal cycle is the 
direct generalization of the unit normal bundle of a smooth submanifold. Its 
formal definition was given in [T4| . Although the normal cycle cannot be defined 
on any compact subset of M^, it exists for a large class of subsets, as convex 
subsets, polyhedra, subsets of positive reach, subanalytic subsets for instance. 
Following mi, any compact subset Q of such that N(0) exists is said to 
be geometric^ and N(0) is called its normal cycle. One of the main properties 
of the normal cycle for our purpose is its additivity [T4|: 

Proposition 1. If Qi and Q 2 are geometric, then Q\ U Q 2 and Qi fl Q 2 are 
geometric and 


N(0i U 02) = N(0i) + N(02) - N(0i n 02). (19) 


Here are some classical examples: 

1. The normal cycle of a smooth submanifold of a Riemannian manifold is 
the closed current associated with its unit normal bundle. 

2. If T) is a compact domain whose boundary is a smooth hypersurface, then 
its normal cycle is the closed current associated with its outward unit 
normal vector field. 

3. If C is a convex body, then its normal cycle is the closed current associated 
with the oriented set 

{(m,e :mGdC,^G ll^ll = e C,< ^,^ >< 0}. 

4. The normal cycle of a polyhedron of can be computed by applying |l9| 
to a decomposition of the polyhedron into (convex) simplices and usingW 


3.2 Asymptotic measures, asymptotic cones 

Let us now define an asymptotic (signed) Radon measure on {resp. TM^) 
associated with a geometric subset. We denote by {resp. Btm^) class 

of Borel subsets of {resp. TM^) with compact closure. 

Definition 2. Let Q be a geometric subset of . 
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• The asymptotic measure defined on TM^, assoeiated with Q and X G 
TL{M) is the map 

I ^ ^ 

defined as follows: 

VBGBTM-,<^g{B)=<-N{g),XBh^ > . ( 21 ) 

• The asymptotic measure defined on , assoeiated with Q and x G 
E{M^) is the map defined as follows: 

yB G =< N(0),x^-i^^^h" > . (22) 

If Q and B are fixed, the map 

is quadratic, inducing its isotropic cone. This remark leads to the following 
definition: 

Definition 3. Let Q be a geometrie subset of . With any Borel subset 
B G , we assoeiate the eone 


eg = {x e S(M^) : = 0} (23) 

and the eone 

= {a; e S(M^) ; a; parallel, ^{3) = 0}, (24) 

whieh we eall the asymptotic cones assoeiated with B. 

In many applications, and for simplicity, it is easier to consider the cone of 
parallel vector fields, identifying a parallel vector field with its value at any point 
m G B. Obviously, C C£. We remark, however, that this new definition 

can be quite restrictive depending on the geometry of . For instance, if 
has non-zero constant sectional curvature, the only parallel vector field is the 
null vector field. In contrast, if = E^, the space of parallel vector fields 
is the space of constant vector fields, which can be identified with E^, which 
is much easier to manipulate. We deduce easily from explicit expressions of 
these curvature measures in some particular cases. 

3.2.1 The case of smooth submanifolds 

Let be a (compact smooth) submanifold (with or without a boundary) 
embedded in and X G TL{TM^). Since the normal cycle of a smooth 
submanifold is its unit normal bundle, we deduce from [8] or [20] (Corollary 16 
page 215) that for any B G Btm^^ 

^^n{B) = / h^ (j)Tj^^nd7rM{X)^pTj^y^nd7rM{X))dfdv^ (25) 

JsT^w^nB 

where h^ denotes the second fundamental form of in the direction of the unit 
vector ST-^W'^ denotes the unit normal bundle of and pr^^^^n denotes 
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the orthogonal projection onto the tangent bundle TW^. In particular, let 
be a (smooth oriented) hypersurface of bounding a domain D. (This 
assumption is not restrictive in our case, since our results are local. It allows 
to simplify some technical points by considering “only one side” of the normal 
cycles (the one corresponding to the outward unit normals), we have, for any 
B in Bmn , 

= / h{pTj^y^N-id7rM{X),pTj^y^N-id7rM{X))dv, (26) 

where ^ is the outward (with respect to D) unit normal vector field of 

and h is the second fundamental form of in the direction We have 

then a correct generalization of|^ Consequently, 

• If 5 is reduced to a point m, 

= {x G TmM^ : hm{prTW^x,prTw^x) = 0}. (27) 

We deduce that is the cone spanned by and that is, 

we have the direct generalization of the corresponding cone defined for 
surfaces in in Section [l] 

• If is a (smooth-oriented) hypersurface of E^ bounding a domain 

then 

= {x e TM^ : / h{prr]^^N-ix,prj^^N-ix)dv = 0}, (28) 

and 

^par,D _ ^ E^ : / h{prr]^^N-ix,prr]^^N-ix)dv = {)}. (29) 

Jw^-^nB 


3.2.2 The case of polyhedra 

to polyhedra. Let P be a domain in E^ bounded by 


We will extend 28 and 


29 


X 


di {N — 1)-dimensional polyhedron For X G 1-L{P^), we can evaluate 

above each simplex. In particular, if is a {N — 2)-simplex, the support 

of N{'D)\^n -2 is the product x Ccr, where C is a portion of circle. Let 

(ei,..., eAr- 2 ) be an orthonormal frame field tangent to Any point of 


^N-2 


X Co- is a couple (m,eAr-i), where m is a point of and e^v-i is a 

unit vector orthogonal to With these notations, we deduce from [8] or 

[20] (Theorem 72 page 216) that for any B G 




{B) = 


E 


1 


-2^QpNJ{cr^-^nB)xC 




We also deduce that for any B e 
= {a; e : / 

J(< 


{a^-^nB)xC 


< X > = 0}. 


(30) 


(31) 


In particular, if I> is a domain of E^ bounded by a polyhedron P^, 31 can be 
reduced to an explicit simple expression: First of all, we identify the (vector) 
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plane e-^ orthogonal to any (oriented) edge e of P with C, as follows : Let 
Til ^ {resp. 712 G e-^) be the unit (oriented) normal to the faces /i, {resp. / 2 ) 
incident to e. Let e+ {resp. e~) be the (oriented) normalized vectors spanning 
the bisectors of ni and 712 (so that (e+,e“,e) is a direct frame of E^). Any 


vector ae+ + be of is now identified with the complex number a 
explicit integration over each term of type (cr^ DP) x C in 


31 


ib. An 
gives the following 


expression of the asymptotic measure and asymptotic cone: 


Proposition 2. 1. For any B G Sjgs and any eonstant veetor field x ofW?, 


eeE 


l{ef\B) 


[{Z{e)\\pr^±x\\‘^ + sm{Z{e))n{{pr^±xf)], (32) 


where E denotes the set of edges of P‘^, Z(e) the angle of the normal to 
the faees ineident to e (being positive if and only if e is eonvex), pr^i. the 
orthogonal projeetion on and IZ^pr^^x^ the real part of the eomplex 
number (prg^x)^. 

2. In partieular, 

gpar.B ^ ^ ]E3 . ^ l{enB) | |2 ^ Q}. 

eGE 

(33) 

We remark that if the Z(e)’s are “small enough” (this can happen for instance 
when P approximates “smoothly” a smooth surface), then sin(Z(e)) is “close 
to” Z(e) and 

$f,(B) - y]^(enB)Z(e) < a:,e- >V (34) 

eeE 

And then, 


- {:c G EV ^(e n B)Z(e) < x, e” > 2 = 0}. (35) 

eeE 


After choosing a scale r, we construct from the previous construction, a eone 
subbundle of TM^ associated with a geometric subset of . Let us denote by 
B{m^r) the ball of radius r, centered at m G . With each point m G M^, 
and for a fixed (small enough) real number r > 0, we associate the cone r)' 

Definition 4. We eall U n the cone subbundle of TE^ at seale r 

meE B{m,r) J 

assoeiated with Q. 


We remark that the dimension of each fiber may change with m. This bundle 
is defined over the whole M^, even at the points m which are “far” from Q. 
If B{m^r) does not intersect the support of then . This 

phenomenon is visualized in Figure]^ where the cone la has its vertex (the black 


point) on the catenoid, and the concern has its vertex out of the catenoid. 
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(a) (b) (c) 



Figure 1: Some asymptotic cones built on a smooth surface (here on a portion 
of a catenoid) (a) An asymptotic cone in green whose vertex is at the center of 
its corresponding Borel set (a transparent ball) (b) The center of the ball may 
be out of the surface (c) When the Borel set is reduced to a single point on the 
surface, the asymptotic cone degenerates to the union of two planes 




Figure 2: Asymptotic cones built on a non-smooth surface, here a triangulation: 
the top of the Lilium tower (designed by architect Zaha Hadid) in Warsaw, 
Poland. 




Figure 3: Asymptotic cones built on a non smooth-surface. Here another tri¬ 
angulation: the Heydar Aliyev Center (designed by architect Zaha Hadid) in 
Baku, Azerbaijan 
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4 Convergence and approximation results 


In the previous paragraphs, we gave our first justification of the denomination 
of asymptotic cones. For surfaces in this cone reduces at each point m 
of to (the product of the normal line by) the standard asymptotic directions 

of at m. We give now a second justification in terms of the convergence of 

sequences of polyhedra of E^. For simplicity, we restrict our study to (oriented) 
smooth hypersurfaces or polyhedra of E^ bounding a (compact) domain. We 
will show, in particular, that, if a sequence of domains Vj^ whose boundaries 
are polyhedra converges to a domain D whose boundary is a smooth 

hypersurface (in a sense that will be clarified later), then for a large 

class of Borel subsets 5, the sequence of asymptotic cones (resp. 
converges to (resp. ). We will use the following terminology m- 




The fatness 0(P^“^) of a polyhedron is defined as follows: If a 

is an /-simplex, we begin to define the size e(cr) of a: It is the maximum 
over all edges e of cr of the length of e. Then, the fatness of a is the real 
number 


©(cr) 


min min 

fj, /-simplex in a 


volj{/a) 

e(cr)J 


Finally, the fatness of P is the minimum of the fatness of its simplices. We 
denote by the class of polyhedra in E^ with fatness greater or equal 
to 0 . 


• An {N — 1) dimensional polyhedron p^-^ in E^ is closely inscribed in 

a smooth hypersurface if its vertices belong to and if the 

orthogonal projection of P^~^ onto is a bijection. 

• Let P^~^ be an {N — 1) dimensional polyhedron (bounding a domain V) 

closely inscribed in a smooth hypersurface in E^. The angular 

deviation am between m G P^~^ and pr^N-i{m) is the maximal angle 
between the normal and {m^Um)^ where {m^Um) 

belongs to the support of N(P). If B is any Borel subset of E^, we write 


aB = sup am- 

mEB 


One of the classical observations is that a “good” fatness of a polyhedron 
and a “small” Hausdorff distance of this polyhedron to a smooth hypersurface 
in which it is closely inscribed imply that the angular deviation is “small”. 


4.1 A convergence result 

Let us state now our convergence theorem. For any cone C, we denote by UC 
the basis of C; that is, UC = C nS^“^(0, 1), where S^“^(0, 1) is the unit sphere 
centered at the origin. We also use the distance d defined on the class of subsets 
of E^ by 

d{A,B)= inf dix^y). (36) 

xeA,yeB 

Theorem 1. Let D he a (compact) domain o/E^ hounded by a smooth hyper¬ 
surface Let {Vk) be a sequence of domains ofK^ hounded by polyhedra 

{P^~^) closely inscribed in such that: 
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1. The limit of {P^ is ^ for the Hausdorff distanee. 

2. The fatness of is uniformly bounded from below by a non-negative 

eonstant: there exists 0 > 0 sueh that for all /c G N; P^~^ G F^. 

Let B G sueh that for all x G = 0. Then, every sequenee 

{xk G of unit veetors admits a subsequenee (still denoted by {xk G C^^)) 

that eonverges to a unit veetor of . In partieular, 

lim = 0. (37) 

k^oo 

In the smooth case, the assumption on the boundary of B can be translated 
in terms of the second fundamental form of W^~^. We say that the normal 
eurvature at a point p of a submanifold V of is null if the second funda¬ 
mental form hp of satisfies hp{u,u) = 0 for every u tangent to V at p. 

We get: 

Corollary 1. Under assumptions^^ and^ of Theorem\^ suppose that B is 
a Borel set sueh that fl dB is smooth and with null normal eurvature. 

Then, the eonelusion of Theorem^ holds. 

To prove this corollary, we simply remark that under these assumptions, 

^^{dB) = / h{prTw^-^^^P'^TW^-^^) = (38) 

Jw^-^ndB 

The rest of this section focuses on the proof of Theorem 

4.1.1 Convergence of sequences of normal cycles 

We need to introduce the flat norm on P/(E^) as follows. The mass of an 
/-current T is the real number 

M(T) = sup{T(c^)}, (39) 

where the supremum is taken over all /-differential forms with compact support 
such that sup^^jgiv \uJm\ < 1. The flat norm of an /-current T is the real number 

jr(T) = inf{M(A) + M(5)}, (40) 

where the infimum is taken over all rectifiable /-currents A and (/ -1- l)-currents 
B such that T = A-\- dB. Our main ingredient in our study of convergence and 
approximation of the asymptotic cones is the following result, which is a simple 
reformulation of Theorem 67 of [20] (page 200) for polyhedra: 

Theorem 2. If P^~^ is a elosed {N — 1) dimensional polyhedron bounding a 
domain D and elosely inseribed in a smooth elosed hypersurfaee o/E^ 

bounding a domain V, then for any Borel subset B of P^~^, 

<A-(4b +afl)M(N(D)|^^j.»), 

where Sb is the Hausdorff distanee between B and pr^N-i{B) and K is a eon¬ 
stant depending on the norm of the seeond fundamental form ofW^~^. 
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The following proposition can be deduced from Theorem in a slightly 
different version, see also m- 

Proposition 3. Under the assumptions^^ and^ of Theorem\^ 

1. The masses M((N(P/c))) are uniformly hounded from above; 

2. The sequenee (N(P/c)) eonverges to N(D) for the flat norm. 

4.1.2 Convergence of sequences of asymptotic measures 

Our framework is the space of (signed) Radon measures on with finite total 
variation, endowed with the norm ||.||i defined for every /i (with finite total 
variation |/i|) by 

Mh = 

It is well known that this space is the (topological) dual to the space Cc(E^) of 
continuous functions with compact support on E^, endowed with the norm ||.|| 
defined by 

11 / 11 = snp \f{x)\. 

The space of (signed) Radon measures (with finite total variation) can also be 
endowed with the topology of weak eonvergenee of measures: a sequence of 
Radon measures (/i/e) on E^ (weakly) converges to /i if, for every continuous 
function / with compact support on E^ {resp. TE^), iik{f) converges to /i(/). 

Proposition 4. Under the assumptions^^ a?id[^ of Theorem\^ 

1. For eaeh veetor x G E^, the sequenee of measures converges to 

for the weak eonvergenee of measures on E^; 

2 . For eaeh unit veetor x G E^, the sequenee of measures (^|)^) is H-Hi- 
hounded. 

Proof of Proposition!^ Item[^is a direct consequence of Proposition!^ since 
the fiat convergence of the sequence of normal cycles (N(P/e)) implies the weak 
convergence of the measures (^x)^)- Item is an application of the Theorem 
of Banach-Steinhaus: Since (^1)^) converges to for the weak convergence 
of measures on E^, for each / G Cc(E^), sup/. | < f >< Too. The 

Theorem of Banach-Steinhaus then implies that 111 uniformly bounded 
with respect to k; that is, the sequence of measures is 11 .| |i-bounded. 

Let us now explain our assumption on the boundary of the Borel subset B 
in TheoremIn general, the weak convergence of the sequence (^x)^) 
imply the convergence of to (4>^(5)) for every Borel subset. Indeed, 

generically, characteristic functions are not continuous. That is why we restrict 
our study to a class of Borel subsets with suitable boundaries with respect to 
such that we can use the following general lemma (see m Chapter 4 for 
instance): 

Lemma 1. Let {fik) he a sequenee of (signed) Radon measures on E^ sueh that 
1. (/i/c) eonverges to fi for the weak topology, 
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2. the sequence (|/i/c|) of total variation of {jUk) converges to a Radon measure 
u, for the weak topology. 

If the boundary dB of B G B^n satisfies u{dB) = 0; then 

Urn = n{B). (42) 

k^oo 


Since the sequence is 11.||i-bounded, we can extract a subsequence 

that converges. From Proposition and Lemma we deduce: 

Proposition 5. Under assumptions^ and of Theorem^ if x e E{E^) and 
B G SjgAT satisfy \^f)\{dB) = 0^ then 

hm = (43) 


The last step of the proof of Theorem is to relate the behavior of the 
sequence for any x to the behavior of their associated quadratic cones. 

For a fixed B G B-^n, we will study the quadratic forms x (B) and 


X ^f){B) introduced in Section 3.2 We use the norm of uniform convergence 


on the space of quadratic forms: A sequence of quadratic forms (pk) defined on 
converges to a quadratic form q if sup||^||^]^ \(lk{x) — q{x)\ tends to 0 when 
k tends to infinity. 


Lemma 2. Let (pk) be a sequence of quadratic forms defined on E^^ which 
converges to a quadratic form q. Let (x^) be a sequence of unit vectors in 
E^^ such that for each k ^ N, Xk belongs to the isotropic cone Ck of pk (I.e. 
Qki^k) = Then there exists a subsequence of (x^) that converges to a unit 
vector X belonging to the isotropic cone C of q. 


Proof of Lemma [21 We have 


lim sup \q{z) - qk{z)\ = 0. (44) 

\\z\\ = i 

Suppose that (xk) is a sequence of unit vectors such that for all k, Xk G Ck- 
Then, by the compacity of the unit sphere, there exists a subsequence (that we 
still denote by (x/^)), that converges to a unit vector x. From[^ we deduce 

lim \q{xk) - qk{xk)\ = hm \q{xk)\ = \q{x)\ = 0, 

k^oo k^oo 

which means that x G C. 

We remark that under the assumption of Lemma we cannot claim that 
the sequence of Hausdorff distances between UCk and UC tends to 0 when k 
tends to infinity, as shown in the following example. Consider the quadratic 
forms in E^ defined for every /c G N* by 

qk{u,v,w) = ^{u‘^ -\-v‘^ Ew‘^). (45) 

We have for every k G N*C/c = {0}, C = E^, UCk = 0 and UC = S^“^(0,1). 

The proof of Theorem follows from Proposition Lemma Proposition 
[Hand Lemma [2 
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4.2 An approximation result 

In this section, for simplicity, we restrict our study to surfaces in We assume 
that (bounding B) and P (bounding V) are fixed, with P being closely 
inscribed in W^. We suppose (without any restriction) that P is endowed with 
a triangulation, denoting by t a generic triangle, and by r{t) its circum-radius. 
The following theorem compares the asymptotic cone of V over a Borel set B 
composed of a union of triangles of P and the asymptotic cone of D over the 
orthogonal projection prw‘^{B) of B on W^. If C is the isotropic cone of a 
quadratic form q, we denote by A^{UC) the set of unit vectors x e-close to (7; 
that is, \q{x)\ < e. 

Theorem 3. Let D he a domain o/E^ hounded hy a smooth hypersurfaee . 
Let V he a domain o/E^ hounded hy a polyhedron P G ^ > 0 elosely inserihed 
in . For any e > 0, there exists rj > 0 sueh that z/max{r(t),t G B} < rj, 
then 

In other words, roughly speaking, under the assumptions of Theorem the 
asymptotic cones of V are close to the asymptotic cones of D. The proof uses the 
results of Section 5-2 of [8]. We summarize them in the following proposition. 

Proposition 6. Under the assumptions of Theorem\^ for every unit veetor x, 

- ^?)(prw^(5))| = Kmax{r(t),t G P}, (46) 


where K a eonstant depending on the area of P, the length of its boundary, 
the geometry ofW^, and 0. 


If X is a unit vector belonging to , then 46 implies that 


|^f)(P)|=Kmax{r(t),tGP}. 


Consequently, if the triangles of P have a sufficiently small circumradius, then 


\^d{vt:w‘^{B))\ < e. 


The conclusion follows. 


5 Some experiments 

5.1 Construction of asymptotic directions of a triangula¬ 
tion 

To mimic the smooth situation, it may be interesting to deduce asymptotic 
directions and asymptotic lines from the asymptotic cones defined on a singular 
surface. This construction may be achieved if one has a natural plane field 
on this surface. This is the case for triangulated surfaces, where each triangle 
spans a plane. Then, to build asymptotic directions at a point of a triangulated 
surface P^, associated to a Borel set P, one can proceed as follows: 

• Consider a point m on P^ (bounding V) and build a Borel set P around 
m; for instance, a ball centered at m with a suitable radius. 
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• Build the asymptotic cone whose vertex is m, associated with B. 

• Build the intersection of with the plane spanned by the face of 

the triangle that contains m. The result is two directions, called the 
asymptotic directions of at m. 

• When m runs over the surface, build a cross field (reduced to a point when 
the cone is reduced to a point). By integrating the cross field, we have 
asymptotic lines. 



(a) (b) 

Figure 4: (a) At each point where the asymptotic cone is not reduced to {0}, we 
build two asymptotic directions on the top of the Lilium tower (b) By integration 
we have asymptotic lines. 



(a) (b) 

Figure 5: (a) The cross field of asymptotic directions on the Heydar Aliyev 
Center in Baku (b) The asymptotic lines obtained by integrating the asymptotic 
directions 

5.2 Approximation of the asymptotic lines of a smooth 
surface 

Using the construction of asymptotic cones, we can approximate the asymptotic 
directions {resp. lines) of a smooth surface in E^. We give an explicit 
example here. 
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First of all, let us consider a portion of a (smooth) catenoid and a 
triangulation closely inscribed in with a sufficiently dense set of vertices. 
As shown in Figures 6a, 6b, 6c, 6d, the intersection of the tangent plane of 


at a point m with the asymptotic cone of a ball centered at m is reduced to 
two lines that are a discrete approximation of the asymptotic directions of 
at m. If necessary, we can approximate the tangent plane itself by the plane 
spanned by a face of the triangulation. By integrating the directions field, we 
obtain discrete asymptotic lines, as shown in Figures and The reader 


can compare the asymptotic lines directly built on the smooth catenoid (these 
lines are orthogonal since the catenoid is a minimal surface) with the “discrete 
ones” obtained by this procedure. 



(a) 





(c) (d) 

Figure 6: (a) Step 1: A green asymptotic cone built at a point m of a triangu¬ 
lated surface approximating a (smooth) catenoid (b) Step 2: The blue plane is 
an approximation of the tangent plane of the (smooth) catenoid at m (c) Step 3: 
The intersection of the green asymptotic cone and the blue plane gives two lines 
intersecting at m (d) Step 4: When m runs over the triangulation, we obtain a 
cross field that approximates the field of asymptotic directions of the (smooth) 
catenoid 
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(a) (b) 

Figure 7: (a) Integrating the cross field gives an approximation of the asymptotic 
lines of the catenoid (b) The cross field is an approximation of the field of 
asymptotic directions on the (smooth) catenoid. They can be compared with 
the blue lines, which are the asymptotic lines directly computed on the (smooth) 
catenoid 

5.3 Comparison of asymptotic lines 

We remark that the previous example gives (roughly speaking) “very good” 
results because the triangulated polyhedron approximating the catenoid suits it 
correctly, in the sense that the “good” fatness of the triangles implies that the 
normal of any triangle t is close to the normals of the orthogonal projection of t 
on the catenoid. To be more precise, we can estimate the error er as follows: Let 
be a smooth surface approximated by a triangulated polyhedron closely 
inscribed on it. Let m be a vertex of P^, and let us denote by ei and 62 the 
asymptotic directions of at m, and by ei, 62 the approximated asymptotic 
directions at m. We define 

eVm = inf (i(Z(ei,ei) + Z(e 2 ,e 2 )), i(Z(ei, 62 ) + Z(e 2 , ei))), (47) 

where all the angles belong to ( 0 , 

For instance, let us consider the portion of the Enneper surface shown in 
Figure and a triangulated polyhedron closely inscribed on it. The error er is 
always less than or equal to 5 degrees. 
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(a) (b) 

Figure 8: (a) Comparison of discrete and smooth asymptotic lines on the En- 
neper surface (b) Some asymptotic lines computed on the polyhedron 

In the following example, we show that, in accordance with the theory, the 
error may be large even with a very thin triangulation closely inscribed in a 
smooth surface, if the angle between the tangent plane at a point and the 
corresponding triangle is too large (the same phenomenon appears when one 
compares the area of a cylinder with the area of a Lantern of Schwarz inscribed 
on it, see [20] for instance). Here, we consider a (smooth) surface in 
obtained as the graph of a tensor product H-spline function 

/ : G = [0,1] X [0,1] ^ M, 

of degree 3 in each variable, defined by 5 x 5 control points over a regular grid 
of 80 X 80 points. Each square of G is triangulated by taking both diagonals. 
We build the corresponding (piecewise linear) triangulation inscribed on 
W^. Then, we compare, at each vertex m of the triangulation, the normal of 
at m and the average of the normals of the triangles incident to m. The 
error varies between 0 and 0.5 degree (see Eigure [91^. On the other hand, we 
compute on the same triangulation the asymptotic directions of and the 
asymptotic cones of P^, from which we deduce discrete asymptotic directions 
by intersecting the cones with the tangent planes as before. Then, we compare 
at each point m the error evm given in Eigurej^ The error at each vertex m 
varies between 0 and 5 degrees, according to the behavior of the normal of the 
faces incident to m with respect to the normal of at m. 
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0.5 


5.0 




Figure 9: (a) Comparison of the normals computed on a B-spline surface and 
an approximating triangulation (b) Comparison of the asymptotic directions by 
computing er^ at each vertex m (the black points are convex points, where the 
asymptotic cone is reduced to {0} (c) Some asymptotic lines computed on the 
triangulation 


5.4 Deformation of asymptotic lines of discrete surfaces 

In the following example, we produce a deformation of “discrete” asymptotic 
lines as follows: We build a triangulation closely inscribed on a smooth surface 
(here a Chen’s surface). We then compute asymptotic lines by the previous 
process using balls of radius R = 3 (the normalization is such that the average 
length of the edges is 1). Then, we slightly modify the position of the vertices 
that can now be out of (in other words, we create noisy data), without 
modifying the topology of the triangulation. With this new triangulation, we 
compute new asymptotic lines using the same process and using balls of the 
same radius {R = 3) or a different radius {R = 6). 
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(a) (b) (c) 

Figure 10: By perturbing the positions of the vertices of the approximated 
Chen’s surface, we get noisy data and corresponding noisy discrete asymptotic 
lines (a) The initial triangulation with asymptotic lines {R = 3) (b) Asymptotic 
lines with noisy data and R = 3 (c) Asymptotic lines with noisy data and R = 6 


6 Conclusion and perspective 

In this article, with any smooth subspace or singular geometric subspace IF of a 
Riemannian manifold M, we associate a family of cones, defined over any Borel 
subset of M. These cones are the generalization of the asymptotic directions 
defined at each point of a smooth surface of the Euclidean space. We obtain 
convergence and approximation theorems, when a sequence of polyhedra tends 
to a smooth subspace. As a consequence, we find good approximations of the 
asymptotic lines of a surface when it is approximated by a suitable triangulation. 
In our future work, we will study the relations between the characteristic of these 
cones and the geometry (and topology) of W. We will apply these results in 
different fields (face recognition to detect similarities between scanned faces for 
instance). 
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